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{*mJ We investigate the newtonian stationary accretion of a polytropic perfect fluid onto a 

, central body with a hard surface. The selfgravitation of the fluid and its interaction with 

luminosity is included in the model. We find that for a given luminosity, asymptotic mass 
q ' and temperature of the fluid there exist two solutions with different cores. 
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1. Introduction 



2. The Shakura model 



The question we want to address in this paper is the following inverse problem: 
having a complete set of data describing a compact body immersed in a spherically 
s ! , symmetric accreting fluid, find the mass of the central body. We assume that we 

■ know the total mass, luminosity, asymptotic temperature, the equation of state of 

. the accreting gas and the gravitational potential at the surface of the core. 

The fundamental question is whether observers can distinguish between gravas- 
tars [1] versus black holes as engines of luminous accreting systems (see a controversy 
in [2, 3]). While we do not address this problem here, we show a related ambiguity 
*-£h ' in a simple newtonian model. 

6 

The first investigation of stationary accretion of spherically symmetric fluids, includ- 
ing luminosity close to the Eddington limit, was provided by Shakura [4]. It was 
later extended to models including the gas pressure, its selfgravity and relativistic 
effects [5-8]. 

In the following we will denote the areal velocity by U(r,t) — dtR (where t is 
comoving time and R the areal radius), the local, Eddington and total luminosities 
by L(R), Le and Lq, quasilocal mass by m{R) and total by M, pressure by p, 
the baryonic mass density by g (the polytropic equation of state will be p = Kg T , 
1 < T < 5/3) and the gravitational potential by <j)(R). The radius of the central 
body is Ro while its "modified radius" is defined by Rq = GM/\</)(Rq)\. Under the 
assumption that at the outer boundary of the fluid the following holds true: 

2 GmjRoo) 2 
Uoo < 5 < floo, (!) 

Woo 

we have the following set of equations: 

M = -AnR 2 gU, (2) 

TTn TT Gm(R) 1_ L(R) , , 
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d R M = 0, (4) 

■ fa 2 a 2 U 2 \ 

L L(R) = M -—--^-^ ftR) j . ( 5 ) 

a is a dimensional constant a — ut / ' 4~Kvn v c. The details of solving the system are 
provided elsewhere [9, 10] and we will present only the main results here. 

We assume that the accretion is critical, i.e., there exists a sonic point, where 
the speed of accreting gas U is equal to the speed of sound a. All values measured 
at that point will be denoted with an asterisk. We define: 

Lq m, Rq 

X= L^ J 7 "#? 
and obtain the total luminosity: 

M 3 ( 2 \ ( 5 ~ 3r )/ 2 ( r - 1 ) 

£o = 0oXooGV^3- (l-y)[y-xeM-xi)} 2 (^fj ■ ( 6 ) 

Xoo is approximately the inverse of the volume of the gas located outside the sonic 
point. In sake of brevity we will use 

iM 2 / o \ (5-3r)/2(r-i) 



(3 = afioXooGTr 2 ^ . . 



to obtain Eq. ([6]) in a form using the relative luminosity: 

x = (3(1 - y)[y - x exp(-x-f)} 2 . (7) 



3. Bifurcation 

For the relative luminosity, fulfilling Eq. |7]) we proved the following theorem: 

(i) For the functional F(x, y) = x — (3(1 — y)[y — xexp(— xj)] 2 there exists a critical 
point x = a, y — b such that F(a, b) = and d y F(a, b) = with < a < b < 1 
and a = 4(3(1 - b) 3 , b= [2 + a exp(-07)]/3. 

(ii) For any < x < a there exist two solutions y(x)~t bifurcating from (a, b). They 
are locally approximated by: 

1/2 

y + = b± f (a -x)[b + aexp(-a7)(l - 2o 7 )] \ 
\(3[b — aexp(— a 7 )][l — acxp(— aj)]J 

(iii) The relative luminosity x is extremizcd at the critical point (a,b). 



4. Discussion 

In the paper we have assumed the existence of an accreting system which satis- 
fies certain conditions. Under those assumptions the complicated set of integro- 
differential nonlinear equations (J2] ESJ) can be simplified to an algebraic one (0). We 
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checked numerically that the performed simplification causes errors of the order of 
1(T 3 (sec [9] for details). 

The analysis of Eq. ([7]) shows that there exist two different solutions, having the 
same total luminosity and total mass, but different masses of the core objects. One 
can also conclude that for sufficiently large (3 the maximal relative luminosity a can 
get close to 1, i.e., the total luminosity approaches the Eddington limit. 

As the two solution branches bifurcate from the point (a, b), there is no much 
difference between the central masses of bright objects (see [9, 10] for plots). How- 
ever, when luminosity is small (Lq <C Le), this difference can become arbitrarily 
large. This can be understood intuitively, because the radiation is small for test 
fluids (since the layer of gas is thin), or when the central object is light (therefore 
weakly attracting surrounding gas). 

The results obtained here are consistent with relativistic analysis neglecting 
interaction between the gas and the radiation [11, 12]. 
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